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Abstract 

Construction of a nonlinear higher-order thermo-hydrodynamics, including correlations, in the 
framework of a Generalized Nonequilibrium Statistical Grand-Canonical Ensemble is presented. 
In that way it is provided a particular formalism for the coupling and simultaneous treatment of 
the kinetics and hydrodynamic levels of description. It is based on a complete thermostatistical 
approach in terms of the densities of energy and of matter and their fluxes of all orders, as well as 
on their direct and cross correlations, covering systems arbitrarily driven- away- from equilibrium. 
The set of coupled nonlinear integro-differential hydrodynamic equations is derived. Illustrations 
of the application of the theory are described in the follow up article. 
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I. INTRODUCTION 



It has been noticed that one of the comphcated problems of the nonequihbrium theory 
of transport processes in dense gases and hquids is the fact that their kinetics and hydro- 
dynamics are intimately coupled and must be treated simultaneously (E.g. see Refs.[l-5]). 
Moreover, on this we may say that microscopic descriptions of hydrodynamics, that is, as- 
sociated to derivation of the kinetic equations from classical or quantum mechanics and 
containing kinetic (transport) coefficients written in terms of correlation functions, is a tra- 
ditional problem of long standing. An important aspect is the derivation of constitutive laws 
which express thermodynamic fluxes (or currents as those of matter and energy) in terms 
of appropriate thermodynamic forces (typically gradients of densities as those of matter 
and energy). In their most general form these laws are nonlocal in space and noninstanta- 
neous in time. The nonlocality is usually dealt with in terms of spatial Fourier transforms, 
and then the laws - now expressed in reciprocal space - become dependent on wave-vector 
Q. A first kinetic-hydro dynamic approach can be considered to be the so-called classical 
(or Onsagerian) thermo-hydrodynamics; it gives foundations to, for example, the classical 
Fourier's and Fick's diffusion laws. But it works under quite restrictive conditions, namely, 
local equilibrium; linear relations between fluxes and thermodynamic forces (meaning weak 
amplitudes in the motion) with Onsager's symmetry laws holding; near homogeneous and 
static movement (meaning that the motion can be well described with basically Fourier com- 
ponents with long wavelengths and low frequencies, and then involves only smooth variation 
in space and time); weak and rapidly regressing fluctuations j^El- 

Hence, more advanced approaches require to lift these restrictions. Consider flrst near 
homogeneity, which implies validity in the limit of long wavelengths (or wavenumber Q ap- 
proaching zero), and to go beyond it is necessary to introduce a proper dependence on Q 
valid, in principle, for intermediate and short wavelengths (intermediate to large wavenum- 
bers). In phenomenological theories this corresponds to go from classical irreversible ther- 
modynamics to extended irreversible thermodynamics 0,1^- This is what has ben called 
generalized hydrodynamics, a question extensively debated for decades by the Statistical Me- 
chanics community. Several approaches have been used, and a description can be consulted 
in Chapter 6 of the classical book on the subject by Boon and Yip Introduction of 

nonlocal effects for describing motions with influence of ever decreasing wavelengths, going 
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towards the very short hmit, has been done in terms of expansions in increasing powers of 
the wavenumber, which consists in what is nowadays sometimes referred to as higher-order 
hydrodynamics (HOH). Attempts to perform such expansions are the so-called Burnett and 
super-Burnett approaches in the case of mass motion, and Guyer-Krumhansl approach in 
the case of propagation of energy (see for example Refs.[ll] and [12]). An usual approach has 
)een based on the moments solution procedure of Boltzmann equation, as in the work of Hess 
isl ] , using a higher-order Chapman- Enskog solution method. The Chapman- Enskog method 
provides a solution to Boltzmann equation consisting of a series in powers of the Knudsen 
number, Kn, given by the ratio between the mean-free path of the particles and the scale of 
change (relevant wavelengths in the motion) of the hydrodynamc fields. Retaining the term 
linear in Kn there follows Navier-Stokes equation, the term in introduces Burnett-like 
contributions, and the higher-order ones {K^ and up) the super-Burnett contributions. 

A satisfactory development of a HOH being also nonlinear and including fluctuations is 
highly desirable for covering a large class of hydrodynamic situations, and, besides its own 
scientific interest, also for obtaining insights into technological-industrial processes having 
an associated economic interest. Also we can mention its fundamental relevance in Oceanog- 
raphy and Metereology, as for example, the study of quite revelant phenomena as thermo- 
haline circulation and ENSO (El Nino Southern Oscillation); see for example Refs.[14] and 
[15] respectively. Moreover, it has been stated [16] that the idea of promoting hydraulics by 
statistical inference is appealing because the complete information about phenomena in hy- 
draulics seldom exists; for example sediment transport, also the more fundamental problem 
in fluid mechanics of describing the velocity distribution in fluids under flow. This latter 
question shall be approached in a future contribution. Indeed, the nonlocal terms become 



specially important in miniaturized devices at submicronic lengths [17[, or in the design of 
stratospheric planes, which fly in rarefled gases in a density regime between the independent 
particle description and the purely continuous description. Another particular problem to 
it related is the one of obtaining the structures of shock waves in fluids for wide ranges of 
Mach numbers . Moreover, Burnett approximation of hydrodynamics has been shown to 
provide substantial improvement on many features of the flow occurring in several problems 
in hydrodynamics, e.g. the case of Poiseuille flow |3| and others |2o| . 

The microscopic derivation of a HOH, together with the analysis of the validity of existing 
theories, is still a point in question. It has been shown ^ that for the case of Maxwellian 
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molecules, whereas Navier-Stockes approximation yields equations which are stable against 
small perturbations, this is not the case when are introduced Burnett contributions to the 
equations. It follows that small perturbations to the solutions, which are periodic in the space 
variable with a wavelength smaller than a critical length, are exponentially unstable. This 
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fact has been called Bohylev's instabilitiy. More recently, Garcia-Colin and collaborators 
have extended Bobylev's analysis for the case of any interaction potential, and have argued 
that one can interpret the fact as to give a bound for a Knudsen number above which the 
Burnett equations are not valid. Moreover, Karlin |23| reconsidered the question looking for 
exact solutions to simplified models: when a linearized ten-moment Grad-method is used, 
and the Chapman-Enskog method is applied to the model, in fact there follow instabilities 
in the higher-order approximations. On the other hand, resorting to the Chapman-Enskog 
solution for linearized Grad ten-moment equations resummed exactly, solutions are obtained 
for which the stability of higher-order hydrodynamics, in various approximations, can be 
discussed. 

Furthemore, inclusion of nonlinearity in the theory, in a Nonlinear Higher-Order Thermo- 
Hydrodynamics (NLHOTH for short and meaning thermal physics of fluid continua), leads 
to additional possible singularities, called hydrodynamic singularities, as, for example, de- 
scribed in Refs.[24] and [25]. A satisfactory construction of a NLHOTH is highly desirable 
for covering a large class of hydrodynamic situations obtaining an understanding of the 
physics involved from the microscopic level, and in the last instance gaining insights into 
technological and industrial processes as in, for instance, hydraulic engineering, food engi- 
neering, soft-matter engineering, etc., which have an associated economic interest. We do 
present here a description of the derivation of a NLHOTH which is based on a formalism 
in Statistical Mechanics shown to be quite appropriate for dealing with systems in far-from- 
equilibrium conditions. 



II. THEORETICAL BACKGROUND 



For building a nonlinear higher-order (generalized) thermo-hydrodynamics on mechanical- 
statistical basis, one needs a nonequilibrium ensemble formalism for open systems. In 
Ref. [26] has been described an information-theoretic approach to the construction of nonequi- 
librium ensembles. It involves a variational method which codifies the derivation of proba- 
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bility distributions - which are also obtained by heuristic approaches or projection operator 
techniques - systematizing the work on the subject of a number of renowned scientists 
pubhshed along the past century. 

According to theory, immediatly after the open system of particles, in contact with 
external sources and reservoirs, has been driven out of equilibrium to describe its state 
requires to introduce all its observables. But this is equivalent to have acess to the so-called 
one-particle (or single-particle), hi, and two-particle, n2, dynamical operators for any subset 
of the particles involved: This is so because all usual observable quantities can be expressed 
at the microscopic mechanical level in terms of these operators (e.g. Refs.[28] and [29]), 
what is described in Appendix A. 

On the basis of the construction of the nonequilibrium statistical operator j^^El, and 
taking into account the fact that a complete description of the nonequilibrium state of 
the system follows from the knowledge of the single- and two-particle density operators 
[or equivalently the density matrix operator, cf. Eqs.(Al) and (A2)], the most complete 
statistical distribution is the one built in terms of the auxiliary ("instantaneously frozen") 
statistical operator 

n{t,0) = p{t,0) X pn, (1) 

where 

p(t,0) = exp|-0(t)-^ j d\ j d'pFix{r,p;t)hix{r,p) 

-Y^jdhj d'p j dh' j c?Vi^2A(r,p,r',p';t)n2A(r,p,r',p')|, (2) 

in the clasical case (the quantum one is given in Appendix B). Index A(= 1, 2, . . . , s) refers to 
the possible different s subsystems of particles (say, different chemical species in a multicom- 
ponent fluid, different classes of quasi-particles in a solid, namely, electrons in Bloch bands, 
phonons, etc.); r and p are the so-called position and momentum field variables. Moreover, 
we have simplified the matter considering the case of contact and interaction of the system 
of interest with ideal reservoir in stationary states, characterized by the statistical operator 
Pr, so that the complete statistical operator can be factorized in the form of Eq.(l). Hence, 
p{t,0) depends on the variables of the system of interest and Pr on the variables of the 
reservoir; both distributions are taken as normalized - as it should -, with 0(t) ensuring the 



normalization of p, meaning that 

(f){t) = y"dr expj -^y^dVy d^pFix{r,p;t)nix{r,p) 

-Y^j d'r j d'p J d'r' J d'p'F^xiv, p, r', p'; t) n^xiv, p, r', p') | , (3) 

and Fix and F2X are the intensive nonequihbrium variables conjugated to ni and n2 (the 
Lagrange multiplcrs in the variational approach). Moreover, dV is the element of the vol- 
ume in the phase space of the system, and for simplicity we have omitted to indicate the 
dependence on F of ni, 77,2, p, 1Z, and that pr depends on the point phase Tr in the phase 
space of the reservoir. 

We recall that p of Eq.(2) is not the statistical operator of the nonequihbrium system, but 
and auxiliary one - as noticed called the "instantaneously frozen quasiequilibrium" statistical 
operator -, but which allows to built the proper nonequihbrium statistical operator (cf. 
Eq.(4) below), which needs to include historicity and irreversibility effects not present in 
p, which does not account for dissipative processes, besides not providing correct average 
values in the calculation of transport coefficients and response functions. 

Finally, the statistical operator explicity written is given by 

n,{t) = exp I In p{t, 0)- J dt' e<''-'^ ^ In p{t', t' - t)^ x pn, 

with p{t, 0) of Eq.(2), and we recall that 

p{t', t' -t)^ exp - t')C^p{t', 0), (5) 

{£ is the Liouvillian operator of the system meaning i£A — {A,H}), which is the aux- 
iliary operator carrying on the mechanical evolution of the system under Hamiltonian 
H ^ Ho + Hi + W ^ Ho + H', where we have introduced H' ^ Hi + W; with ^0 be- 
ing the kinetic energy operator. Hi contains the internal interactions and W accounts for 
the interaction of the system with reservoirs and sources. Finally, e is an infinitesinal positive 
real number which is taken going to zero after the traces in the calculation of averages have 
been performed (it is present in a kernel that introduces irreversibility in the calculations, 
in a Krylov-Bogoliubov sence). We stress that the second contribution in the exponent in 
Eq.(4) accounts for historicity and irreversible behavior from the initial time (taken in the 
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remote past, to —oo, implying in adiabatic coupling of correlations, (see for example 
Ref. [26]), or alternatively, can be seen as the adiabatic coupling of the interactions respon- 
sable for relaxation processes 30|). Moreover we notice that the time derivative in Eq.(4) 
takes care of the change in time of the thermodynamic state of the system (in the first term 
on the argument i.e. t') and of the microscopic mechanical evolution [second term in the 
argument, i.e. t' — t and see Eq.(6)], and that the initial value condition is TZeito) = p(to, 0) 
for to — oo. 

But, as shown elsewhere [26 1 it is quite convenient, and intuitively more satisfactory, to 
work in an alternative description than the one in Eq.(2), namely, the so-called generalized 
nonequilihrium grand- canonical ensemble. In words, it consists into introducing as basic 
variables independent linear combinations of ni and n2. This can be done along several 
ways 27[, but once we are here working within the framework of classical mechanics, we 
present a simple derivation in Appendix C. According to the results presented in the latter, 
and once for simplicity we are considering the case of the presence of only one kind of particles 
in the system [i.e. A = 1 in Eq.(2)], the auxiliary generalized classical nonequilibrium grand- 
canonical statistical operator is given by 

p(t,0) = exp{-5(t,0)}, (6) 

where S{t, 0) is the so-called informational entropy operator S^l, which we write as composed 
of two contributions, namely 



S{t,0) = Si{t,0) + S2{t,0), 



(7) 



with 5*1 and S2 given by 



^i(t,0) = 0i(t)+ / rfV<^ A(r,t)n(r)+F,(r,t)/i(r) + V(r,t)-t(r) + F,(r,t)-i,(r) 



>2 



(8) 



Ti n' 1-/ 



(9) 



p,p' r,r' 



where p and p' stands for indexes n or h, and r and r' are equal to zero (the densities), 1 
(the vectorial fluxes), and 2, 3, . . . (higher order fluxes). In the last two equations we have 



introduce the definitions 



h{r)= I d3pi-ni(r,p), (10) 



n(r) = J d^pni{r,p), (11) 

i,(r) = J rf3p|!-uni(r,p), (12) 

i„(r) = J d^pnni{r,p), (13) 

7M(r) = I d'p^u^r^ Mr,p), (14) 

/M(r) = y d^pKM^^(r,p), (15) 

where [cf. Eqs.(C8) and (C9)] 

li''"' = [u . . . r — times . . . u], (16) 

with [. . .] being the tensorial product of r times the generating velocity (a vector) 

u = p/m, (17) 

which is the group velocity of the particle with energy- dispersion realtion e{p) = p^/2m, 
but the present formalism also encompasses the case of systems with any kind of energy- 
dispersion relation, as, for example, electrons in quantum states |k > of Bloch bands (when 
we take p = ^ k) , or phonons associated to lattice vibrations (for acoustic phonons in Debye 
model £{p) = s |q| (taking p = with q running in Brillouin zone, and where s is the 
velocity of propagation), and similarly for photons in black-body radiation with ^(p) = c |k| 
(for p = hk, and c is the velocity of light): the general case is briefly summarized in 
Appendix D. 

The operators in Eqs.(lO) and (11) are, respectively, the density of kinetic energy and 
of particles, whose integration in space provides the kinetic energy Hamiltonian and the 
number of particles. The vector operators of Eqs.(12) and (13) are the fluxes of the two 



previous densities, called first-order fiuxes or currents, and those of Eqs.(14) and (15), with 
r > 2, are the higher-order (tensorial) fluxes of kinetic energy and of matter (or particles). 



where, we recall, [. . .] stands for the tensorial product, p and p' for h or n (referring to energy 

and particle densities), and r and r' are 0, 1, 2, 

The average values of the operators of Eqs.(lO) to (15) - associated to the single-particle 
dynamical operator -, and those of Eq.(18) - associated to the two-particle dynamical oper- 
ator -, define the set of basic macrovariables for the thermo-hydrodynamics, namely 



for r = (the densities) , r — 1 (the vectorial fiuxes or currents) , and r > 2 (the higher-order 
fiuxes). Hence, we do have 



All of them are present in Si, while in §2 we find the tensorial operators (of rank r -|- r') 



(18) 




(19) 




(20) 




(21) 




(22) 




(23) 




(24) 




(25) 



where 




(26) 
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Similarly for the contributions associated to the two-particle dynamical operator we do 
have the corresponding thermodynamic variables 

cir\vy.t)^ j dT&;^r\vy)n,{ti (27) 

with C given by Eq.(18). 

The complete set of thermo-hydrodynamic variables are presented in Table I organized, 
as one is going up in both columns, in increasing order of tensorial rank. 



Generalized Grand-Canonical Description 


Single-Particle 


Two-Particle 






7M(r,t);/r(r,i) 


Cl;r'\v,v',t)-n + r2 = r 






/f(r,i);llkr,i) 


C(r,r',*);CPri(r,r',t);C£t^](r,r',0 


lh{v,t) ; In(r,t) 


C(r,r',t);C£rl(r,r',t); 


h{Y,t) ; n(r,t) 


C(r,r',*) 


Dynamical Density Operator 


Dynamical Density Operator 


ni(r,p) 


n2(r,p,r',p') 



TABLE I: The basic set of macrovariables one- and two-particle dynamical operators description 
and the generalized grand-canonical description. 



We call the attention to the fact that the generalized nonequilibrium grand- canonical 
statistical operator is the one of Eq.(4) once in it is introduced the auxiliary operator of 
Eq.(7); explicitly 

n,{t) = p,{t) X PR, (28) 

with 

p,(t)=exp|-%0) + J' dt'e<''-'^^S{t,t'-t)y (29) 

Purthemore, we recall that it can be introduced the separation of pe into two parts [coming 
from each contribution in the exponent of Eq.(20)], namely 

p,(t)^p{t,0) + p'M, (30) 
10 



with p of Eq.(6), and that for the basic variables (in this case those of Table I) and only the 
basic variables, it follows that 

I^\v,t) = I dT il^\T) p^it) = I dTil^\T)pit,0), (31) 

Ci;r'kr,r';t) = J dT Clr'\r,r') p,{t) = J dT ClV\,r') -p{t,0), (32) 

i.e., p[ does not contribute for the average value of, we stress, the basic variables only 
Finally, it can be noticed that the quantities of Eq.(32) can be related to fluctuations of the 
densities and their fluxes. In fact we can write 

^!;r'kr,r',t) = JdV - I^;\r,t)) {f;\r')-I^,\r',t)) p,(t) 

= c!;;'\r,r'-t)-4^\r,t)iP{r',t), (33) 

which are instantaneous at time t but containing space correlations. 

Let us next proceed to derive the equations of evolution of the basic variables, that is, 
the generalized hydrodynamic equations. 



III. NONLINEAR HIGHER-ORDER THERMO-HYDRODYNAMICS (NLHOTH) 

We begin considering a NLHOTH for a system of independent particles, say, an ideal 
fluid or, more generally, a system of particles with the interaction between them treated in 
an average-field approximation (a quite interesting case is the fluid of mobile electrons in 
crystals). Hence, the set of basic variables consists of the one on the left side of Table I, and 
Eq.(19), that is the one composed of 

{hir, t), n(r, t), I,(r, t), I„(r, t), {/^ (r, t) } , {/^(r, t) } } . (34) 

The hydrodynamical motion of this fluid of single particles is described by the set of coupled 
highly-nonlinear integro-differential equations consisting of the kinetic equations provided 
by the formalism (see Appendix E), which are the generalized hydrodynamic equations given 
by 

|-/i(r,t) + V-I;,(r,t) = 4(r,t), (35) 
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d_ 
di 



n{r,t) + V ■In{r,t) = Jn{T,t), 



(36) 



d_ 
di 



Uv,t) + V-ll^Kr,t) 



Jft(r,i), 



(37) 



d_ 
di 



I„(r,t) + V-/|Kr,t) 



(38) 



^(r,t) + V-/ri(r,t) 



JVM 



(39) 



d_ 
di 



^r(r,^) + V-/r'l(r,i) 



(40) 



with r > 2, V- is the tensorial divergence operator, and where, on the right, are present the 
coUision operators 



for r = 0, 1, 2, . . ., p = /i or n, {. . .} is Poisson parentheses (see Appendix E), and subindex 
nought stands for evolution in the interaction representation. We recall that H' contains the 
internal interactions plus the one with the reservoir, and, eventually, with external sources of 
perturbation. Moreover, if the collision operators are neglected we do have the conservation 
equation for each quantity. As a matter of practicality, one needs - depending on the problem 
in hands - to resort to a truncation in the chain of generalized hydrodynamic equations, 
Eqs.(35) to (40), say for r = n: the point is discussed in Ref. [26], and further considerations 
are presented in Refs.[ll,27]. Hence, once the equation for the n-th flux contains a term with 
the divergence of the (n+l)-th flux, to close the system of hydrodynamic variables one needs 
to express the latter in terms of all the other basic variables (the densities and the fluxes up 
to order n), and this involves a nonlinear expression. Thus we do face nonlinearities at this 
point, but also, and more importantly, in the collision integrals J. They are - see Eq.(41) 
- depending on a highly nonlinear way on the nonequilibrium thermodynamic variables F's 
[cf. Eqs.(8) and (9)], and therefore the system of hydrodynamic equations needs be coupled 
with the set of nonequilibrium equations of states, namely, Eqs.(20) to (25) and (27) (let 
us noticed that such dependence comes from the statistical operator in the calculation of 
averages). Consequently, the collision operators depend on a highly nonlinear way on the 
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nonequilibrium thermodynamic variables, which, in turn depend on a highly nonlinear way 
on the basic hydrodynamic variables. All this is illustrated in the follow up article in a 
simple model and a first-order hydrodynamics (n = 1). 

Summarizing, Eqs.(35) to (40) provide the thermo-hydrodynamics of the nonequilibrium 
many-body system, which describe the evolution of an infinite set of coupled highly-nonlinear 
integro-differential equations. Moreover, as noticed, practical usage requires us to introduce 
a truncation procedure, meaning the analog of the one in the Hilbert-Chapmann-Enskog ap- 
proach to the Boltzmann equation, or, more precisely, one which is closely related to Grad's 
moments method j^^]: This is discussed in Ref.[33], and a particular case - a phoinjected 
plasma in semiconductors - is considered in Refs.[34] and [35]; we return to this important 
question in a future article on NLHOTH. 

It can also be noticed that, alternatively, Eqs.(35) to (40) can be closed writting on 
the right-hand side the basic variables in terms of the intensive ones, but for building this 
NLHOTH it is more satisfactory the first choice of description (i.e. in term of the variables 
of Eq.(34)). 

We also call the attention to the fact that if we introduce a nonequilibrium generalized 
grand-canonical-like partition function as 

0(t) = lnZ(t), (42) 

with (j){t) of Eq.(3) but for F2 = once we are working in a single-particle representation, it 
follows that 

ll^\r,t) = -6\nZ{t)/6Fp{r,t), (43) 

an alternative of Eq.(23) and where 6 stands for functional derivative j^^. Moreover, defining 
the nonequilibrium-informational entropy 



S{t) = J dTS{t,0), (44) 

it is verified that 

F;,%,t) = 6S{t)/6ll^\r,t). (45) 

It can be noticed that Eq.(45) and Eq.(43) - the latter an alternative form of the Eqs.(20) 
to (25) - relate the basic variables (the densities and their fluxes) Ip\r,t), to the nonequi- 
librium thermodynamic variables, Fp\r,t), which, as noticed, can be labelled, in analogy 
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to the case of systems in equilibrium, nonequilibrium equations of state. As also commented 
above, they are fundamental for closing the generalized hydrodynamic equations, Eqs.(35) 
to (40), that is, the set of equations of evolution is coupled to the set of nonequilibrium 

equations of state. 

On the other hand, if the physical situation we are considering requires to further intro- 
duce as basic variables the pair-correlation-like contributions of Eq.(27), use must be made 
of the statistical operator of Eq.(4) where the auxiliary statistical operator of Eqs.(6) to (9) 
must be introduced. The set of Eqs.(35) to (40) - in which it must be understand that the 
average values are calculated using the above said statistical operator -, is now extended to 
incorporate the equations of evolution of the added variables, namely 



d_ 
di 



Ci;r'kr,r';t) = J dF {clr\r,v'), Ho} pit,0) 

+ jdTj dVn {Ci;r'l(r,r'),^'} ^ x (46) 



or 



d_ 
di 



C(r>r';t)+ / dV [v-i^'+il(r)/;;i(r')+/;l(r)V'-/;:+^kr)] p(t,0) 
= |cj;r'l (r, r'; t) + V • / dV (r, r'; t) + V ■ j dV Cl^^'-^^^ (r, r'; t) 

- jlT'\r,r';t), (47) 

The set of equations of evolution needs, as noticed, be coupled to the equations of state, 
that is, those of Eqs.(43) or (45), and now with the ones associated to the two-particle 
contributions, namely 

C;;r'kr,r';t) = 51nZ(t)/5F[;r'l(r,r';t), (48) 



or 

SF^r\r,r';t) = -SS{t)/SClr\r,r';t). (49) 

We restate that both cases, the one involving only single-particle operators and the other 
with single- and two-particle operators, the nonlinearity is present in the scattering oper- 
ators i7's, which are in general depending on a highly nonlinear variables F's, and then, 
through the nonequilibrium equations of state, depend on a highly nonlinear way on the 
basic variables. This is illustrated in the follow-up article. 
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IV. CONCLUDING REMARKS 



Briefly summarizing the results, it has been shown how a statistical nonequilibrium en- 
semble formalism (applicable to the study of systems even in conditions far from equilib- 
rium) provides - via the use of a quite general nonequilibrium grand-canonical ensemble - 
a microscopic foundation for a Nonlinear Higher- Order and Fluctuations-dependent Hydro- 
dynamics. Its description is based on the set of macrovariables consisting of the densities of 
energy and matter (particles), their fluxes of all order, and the direct and cross correlations 
between all of them: In that way higher-order descriptions and fluctuations, respectively, 
are included in the theory, (cf . Table I) . 

Al these macrovariables are the average value over the nonequilibrium ensemble of the 
corresponding microscopic mechanical operators. Once the complete set of macrovariables 
is given we can obtain the nonlinear hydrodynamic equations, which are the average value 
over the nonequilibrium enesemble of Hamilton equations of motion (in the classical level or 
Heisenberg equations at the quantum level) of the basic microvariables (mechanical observ- 
ables): cf. Eqs.(35) to (40) and (46) to (47). 

These generalized hydrodynamic equations present on the left side the conserving part 
of the corresponding quantity, and on the right-hand side are present the collision integrals 
which include the action of external sources and the contributions of scattering processes 
responsible for dissipative effects. As noticed before, these collision integrals depend on a 
highly nonlinear way on the nonequilibrium thermodynamic variables (the Lagrange mul- 
tipliers that the variational method introduces, namely, F^^' and -Fp^^^^ in Eqs.(7) to (9)). 
Hence, this set of hydrodynamic equations is closed after it is coupled with the set of nonequi- 
librium equations of state, that is, Eqs.(43) and (48), or Eqs.(45) and (49) or simply Eqs.(19) 
and (27). 

In that way we do have a quite generalized hydrodynamics under any arbitrary con- 
dition of excitation, which, as noticed, can be referred to as Nonlinear Higher-Order and 
correlation-dependent Hydrodynamics. In the follow up article we present an illustration of 
its application, consisting of a simple model of two ideal fluids in mutual interaction - one 
of them acting as a reservoir - and in a truncated NLHOH of order 1. 
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APPENDIX A: DYNAMICAL DENSITY OPERATORS (REDUCED DISTRI- 
BUTION FUNCTIONS) 

For the sake of completeness, we notice that in classical mechanics the one-particle and 
two-particle operators fii and ri2 are given, respectively, by 

N 

^i(r,p) = Xl'^('^~'^j)'^(P~Pj)' (^1) 

N N 

Mr, P, r', P') = E E '^(^ - ^^■) '^(P - P^) ^(^' - ^'^) '^(P' - P^)' (A2) 
j=i k=i 

for each one of the different S types of particles that can be present [cf. Eq.(2)], and where, 
as usual, rj and pj are the coordinate and momentum of the j-th particle, and r and p the 
so-called coordinate and momentum field variables. 

In quantum mechanics the one- and two-particle density operators are {a is the spin 
index) 

fii(r,a,r',a') = *t(r), (A3) 

ri2(ri,(7i,r2,(72,r^,a^,r;,(j;) = *t^(ri) *^^(r2) *4(r^)*^i(r;), (A4) 

where ^(^"f") arc single-particle field operators in second quantization (an excellent didactical 
description of them is available in the article by B. Robertson of Rcf.[38]). These operators 
can be expressed in terms of annihilation and creation operators in single-particle states 
with, say, wavefunction (/^kg. and energy levels e^^, in the form 

*^(r) = ^(^kaCka, (A5) 

k 

*Ur) = E<^k.cL, (A6) 

k 
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This is quite convenient since it introduces a description in the reciprocal k-space, with 
practical advantages in calculations and physical interpretations. In this way we can intro- 
duce instead of hi and h2 of Eqs.(A3) and (A4), matrices composed with the quantities 



Ck<7 Ck'fT' 



(A7) 



Hence, any one-particle mechanical observable can be expressed in terms of the quantities 



(A8) 



with Q 7^ 0, where we have separated out diagonal terms and non-diagonal ones: the first 
are so-called populations and the others coherences 123] (while the former appear in the 
description of whole properties - or uniform states -, the latter describe local properties - or 
inhomogeneous states). Any one-particle mechanical observable is then written in second 
quantization in the known form 



A= J2 (ka|i|kV^ 4^Ckv', 

kcr.k'cr' 



(A9) 



for any pair kcr and kV, and where is present the matrix element of A between one-particle 
states ^kcr| and IkV^ 

For two-particle mechanical observables, using in Eq.(A7) the change of indexes k -|- 
for k, k' - iQ' for k', k - iQ for ki, k' + |Q' for kg, we have 



hkai,k'a[,k' ,a2,rk,a'^ — cj^o-^ cj.,^^ Ck'o-2 ^kcr^; 

for Q, Q' 7^ 0. A two-particle observable takes the form 0, llll 



(A 10) 



k,k',Q,Q' cricr[cr2Cr'2 



k + ^Qai, k' - ^QV2|5|k' + ^Q'a'„ k - ^Qa[^ x 



where we have rearranged the order in the product of four operators in Eq.(All), to obtain a 
contribution as a product of two ni-type operators and amother involving hi already present 
in Eq.(A8); upper sign minus stands for the case of fermions and lower plus for bosons. 
For example, Coulomb interaction between pairs of electrons is given by 



hQ hl^ - 



(A12) 
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where, once a complete set of plane-wave states has been taken as one-particle quantum 
states, V(Q) = Ai^e^ /VQ"^ (V is the sample volume), 

= ^ ^aa' ^kQo-a', (Al3) 
ko-o-' 

is the Q-wave-number Fourier transform of the operator density of charge of the electrons 
(in units of the electron charge), and N — X^ko- ^ko- ^ko- is the operator number of particles. 

We call attention to the fact that if the particles are bosons, then coherent states are 
possible to be present and then, besides the one- and two-particle density operators, we 
must also add the individuals amplitudes Cko- (and c\^^) for the description of the boson 
system (for example coherent photons in a laser, coherent phonons in semiconductors, etc). 

APPENDIX B: THE QUANTUM MECHANICAL NONEQUILIBRIUM STATIS- 
TICAL OPERATOR 

In the quantum case the auxiliary statistical operator of Eq.(2) is 

p(t, 0) = exp \ -<P{t) - 5^ i^k(t) nk - J] 5^ FkQ,k'Q' (t) n^Q 4'Q' [ , (Bl) 

Ik kQ k'Q' J 

which is written in terms of the quantities of Eqs.(A8) and (All), with spin indexes omitted. 
The quantum nonequilibrium statistical operator is given by Eq.(4), where now enters p of 
Eq.(6), and 

p{t, t'-t)^ exp |-l(i' - i) ^1 pit, 0) exp |T^(t' - i) ^1 . (B2) 

The three F's in Eq(Bl) are the corresponding Lagrange multipliers in the variational ap- 
proach (intensive nonequilibrium thermodynamic variables conjugated to the dynamical 
operators). Moreover, the macrovariables for describing the nonequilibrium thermodynamic 
state of the system are the average value of the same quantities over the nonequilibrium 
ensemble, namely 

r^k(^) = ^r{r^k7^,(^)}, (B3) 
nkQ(t) = Tr {nkQ7^,(^)}, (B4) 
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with Q 7^ 0, and 

rikQM'Q'it) = Tr |nkQni/Q' T^e{t)^ , (B5) 
with no restriction on Q and Q'. 

APPENDIX C: CLASSICAL NONEQUILIBRIUM GRAND-CANONICAL STA- 
TISTICAL OPERATOR 

Consider thermodynamic intensive variable Fi{r, p; t) in the auxihary statistical operator 
of Eq.(2), and let us perform a Taylor series expansion in p, i.e. 

F^{r,p■,t) = $(r,i) + $(r,i).p + $Pl(r,i)®[pp] 

+ ^ $M(r, t) (g) [p . . . r - times . . . p], (CI) 

where $ = Fi(r, p; t)|p=o, = ^[Vp ... r - times ... Vp]Fi(r, p; t)|p=o, for r = 1,2,3,... 
indicating the tensorial rank, [. . .] stands for tensorial product of vectors, dot as usual as 
scalar product, and <S> for fully contracted product of tensors. Next, we redefine quantities 
$ in the following way 

^r,t)=A^{r,t), (C2) 
^(r,t) = Vn(r,t), (C3) 

$M(r, t) ^ + \ [($M(r, t) ® 1^1) 1^1] , (C4) 

further introducing 

"^V^ = -^^'K^^t) (8) iP], (C5) 
where 1^^^ is the unit tensor of rank two, it can be noticed that it is satisfied the rule 

FM(r,t)® iP] =0, (C6) 

and that 

P/iiM (r, i) [p . . . r - times . . . p] = (r, t) (p) + Pj^'^^ (r, t) ® , (C7) 
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where 

■ut'"' = [u . . . r — times • • • u], (C8) 
is the r-rank tensor resuhing from the inner tensorial product of r characteristics velocities 

u = p/m. (C9) 
Finally, defining the generalized fluxes operators 

/M(r)= / rf=^pMMni(r,p), (CIO) 



/M(r) = I d3pZ.^W^^(r,p), (Cll) 

and taking into account that 

n2(r, p, r', p') = h,{r, p) n,{r', p'), (C12) 
we can write the auxiliary statistical operator of Eq.(6). 

APPENDIX D: THERMO-HYDRO DYNAMIC VARIABLES FOR SYSTEM 
WITH ARBITRARY ENERGY-DISPERSION RELATION 

As shown elsewhere jsgl at a quantum-mechanical description the fluxes of energy (index 
h) and of matter (index n) are in reciprocal space (of wavevectors Q) given by 

^'(Q,^) = E^(^k+iQ + ^k-iQ) «''■'(k,Q)r^kQ(t), (Dl) 

k 

/M(Q,t) = ^nM(k,Q)nkQ(t), (D2) 

k 

where 

uM(k,Q) = [u(k,Q)...r- times... u(k,Q)], (D3) 

is a rank-r tensor composed by the tensorial product of r-times the generating velocity 

1 ~ 1 /I 

u(k.Q) = -V... + |:p^(-Q.V.j -v.... (D4) 
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We recall that is the energy dispersion relation of the single particles, and 

^kQ(t) = Tr Ck_iQ Pit, 0)} , (D5) 

is the average over the nonequilibrium ensemble of Dirac-Landau-Wiegner single-particle 
dynamical operator. It can be noticed in Eq.(D4) that the generating velocity is composed 
of a first contribution consisting of the group velocity of the wavepacket of a single particle 
in state |k), and the following terms - depending on Q - can be considered as contributions 
indicating the change from point to point of the group velocity, once in direct space we can 
write 

u(k, r) = I^Vk^k + ^ (V ■ Vk)' ^Vk^k + . . .} 5(r), (D6) 

where V is the gradient in r-space, and we can see that this expression contains a series of 
even powers of such gradient. These expressions in a classical mechanical levels retain the 
form except that we must introduce the linear momentum p in place of h k. 



APPENDIX E: THE KINETIC THEORY 



The kinetic theory that can be derived out of the nonequilibrium ensemble formalism 
can be considered a fra-reaching generalization of Mori's approach, involving nonlinearity, 
all order of particle collisions, and covering arbitrarily fra-from-equilibrium situations. The 
presentation unavoidably involves a lengthy and somehow cumbersome mathematical han- 
dling, but we call the attention to the simplest form (the Markovian approximation) which 



is usually t 
validate 



le one to be used, once in most problems are verified the restrictions that, say. 



it 3, 42, 3- 

We begin recalling the separation of the total Hamiltonian as given by 

H = Ho + H', (El) 

where H' = Hi + W, with Hi containing the internal interactions and W those with the 
reservoir and eventual external sources of pertubation, while Hq is the kinetic energy oper- 
ator. Moreover, let us call in a generic way as J = I5 2, . . . the set of basic variables, 
which for the case we are considering on the main text are those of Eq.(26), for which stands 
the selection rule 

P,,Ho] =J2(^^rnPm, (E2) 

m 
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where c-numbers or differential operators depending on the represdentation being 

used, and in the case of the variables of Eq.(26) we have that 



(E3) 



with p — h or n and r — 0,1,2, We also recall that for the nonequilibrium statistical 

operator of the system holds the separation into two parts, as given by Eq.(22) namely 



(E4) 



consisting of the "instantaneously frozen" p, and the part p'^ responsible for the irreversible 
evolution of the macrostate of the system. 

The equations of evolution for the basic variables are the average over the nonequilibrium 
ensemble of Hamilton equation of motion for the basic dynamical variables, namely 



d 
di 



Qjit) = jTr {PjPeit) xpR^^Tri^^ [p„ h] p,{t) x p«| , 



(E5) 



where Tr is to be understood as integration over the composite phase space of system and 
reservoir. Using Eqs.(Dl) to (D4) we can alternatively write that 



Pj,Ho 



p{t,Q)xpR}+Tr\ — 



P'eit) X PR 



P„H' p{t,0)xpR\+Tr\- Pj,H' p',{t)xpR 



(E6) 



But the second term on the right-hand side of this Eq. (D6) is null; in fact 



Pj,Ho 

once we take into account that 



Tr I [Pj, Ho] p'eit) xpr\ = J2 Oijm Tr [Pmp'e{t) X Pr} = 0, (E7) 



Tr {P^ pe{t) X pn'j = Tr {p^ p{t, 0) x pn'j , 



(E8) 



[cf. Eq.(23)]. The remaining term can be handled through an elaborate algebra to ob- 
tain that the kinetic equations acquire the form of a far-reaching generalization of Mori's 
equations, namely 



d 
di 



(E9) 



m=2 
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where J*^°^ is the first contribution on the rightof Eq.(D6), J^^^ the third followed by a 
series of partial collision operators involving increasing orders in the interaction strengths 
(m = 1, 2, . . .), and containing memory effects, space correlations, and are highly nonlinear in 
the basic variables (the collision operators on the right-hand side of Eq. (E9) are dependent on 
the intensive nonequilibrium thermodynamic variables, F's in Eq.(2), which are dependent 
on the basic variables through the constitutive equations, cf. Eq.(19) or Eqs.(43) and (45)). 

In the extremely complicated set of coupled equations (E9), we can single out the sim- 
plified form consisting of the markovian approximation. The Markovian limit of the kinetic 
theory is of particular relevance as a result that, for a large class of problems, the interactions 
involved in H' are weak and the use of lowest order in the equations of motion constitutes 
an excellent approximation of good practical value. In chapter 6 in Ref. [26] are described 
several examples of its application for which there follows an excellent agreement between 
the calculation and the experimental data. By means of a different approach, E. B. Davies 
has shown that in fact the Markovian approach can be validated in the weak coupling 
(in the interaction) limit (retaining only the quadratic contribution). 

Making it explicit, the Markovian equations in the kinetic theory are 



d 



dt 



(ElO) 



where we have introduced explicitly the eventual dependence on position r. Th e three terms 
on the right-hand side of Eq.(DlO) are 



Tr < 


in 




, Ho 


Tr < 


Jn 







p(t,0) X PR 



(Ell) 



(E12) 



.^^^(M)^(i- 



2 rt 



dt'e<''-'^ xTr 



H'{t'-t)o, 



pit,0)xpn\, (E13) 



iijfir,t) = ^J2j^ cit'e^(*'-*) xTr[[H'{t'-t)o,Pk{r)]p{t,0)xpn} 



6Qk{r,t) 



and we recall that J-°^ and J]''' , which in Mori's terminology 



r(i) 



are called precession terms. 



are related to the non-dissipative part of the motion, while dissipative effects are accounted 
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for J(^). In many cases, as a result of particular symmetries in H' and p, there follows that 



Finally, we notice that in Eqs.(35) to (40), we do have in the right-hand side that the 
collision operator is 



Jj'^](r,t) = J^'^]«(r,^) + ^^]M(-)(r,t), 



(E15) 



m>2 



in the general case, or 



jJ'^l(r,t)^4^1('Hr,t) + 4'-l(2)(r,t), 



(E16) 



in the Markovian regime j26 
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